Polarization entangled state measurement on a chip 
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The emerging strategy to overcome the limitations of bulk quantum optics consists of taking ad- 
vantage of the robustness and compactness achievable by the integrated waveguide technology. Here 
we report the realization of a directional coupler, fabricated by femtosecond laser waveguide writing, 
acting as an integrated beam splitter able to support polarization encoded qubits. This maskless 
and single step technique allows to realize circular transverse waveguide profiles able to support the 
propagation of Gaussian modes with any polarization state. Using this device, we demonstrate the 
quantum interference with polarization entangled states and singlet state projection. 



Photons are the natural candidate for Quantum In- 
formation (QI) transmission [H, Q, quantum computing 
@, 0] optical quantum sensing and metrology [5] . How- 
ever, the current optical technology does not allow the 
transition to ultimate applications for many practical 
limitations. Complex quantum optical schemes, realized 
in bulk optics suffers from severe limitations, as far as 
stability, precision and physical size are concerned. In- 
deed, it is a difficult task to build advanced interfero- 
metric structure using bulk-optical components with the 
stability and the optical phase control accuracy necessary 
to reach the sensitivity allowed by quantum mechanics. 
Furthermore, it is very difficult to reach this task out- 
side environments with controlled temperature and vi- 
brations, and this makes applications outside laboratory 
hard to achieve. 

The present approach to beat these limitations is to 
adopt miniaturized optical waveguide devices. Very re- 
cently it was reported 0, 0] that silica waveguide cir- 
cuits integrated onto silicon chips can be successfully used 
to realize key components of quantum photonic devices. 
Inherently stable interferometers were shown to demon- 
strate phase stability, not only of single path encoded 
qubits, but also of two-photon entangled Fock state. On 
this basis, miniaturized integrated quantum circuits were 
realized to implement the first integrated linear optical 
control-NOT gate, achieving a fidelity very close to the 
theoretical value [6[ . More recently, novel components for 
adaptive quantum circuits have also been demonstrated 
[i[ . These experiments demonstrate robust and accurate 
phase control in integrated, path-encoded waveguide sys- 
tems. Similar results have been obtained in UV laser 
written optical circuits fabricated in a suitable stack of 
doped silica layers on a silicon substrate [9[. 

All the experiments performed so far with integrated 
quantum circuits are based only on path encoded qubits 
with a given polarization state of the photons. On the 
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FIG. 1: Schematic of the femtosecond- laser- written direc- 
tional coupler in the bulk of a borosilicate glass. Upper in- 
set shows a microscope image of the two waveguides in the 
coupling region. Lower inset shows the near- field intensity 
profile of the output guided modes of the directional coupler 
by launching light in a single input; the symmetric Gaussian 
shape and the balanced splitting in the two arms can be ap- 
preciated. 



other hand, many QI processes and sources of entangled 
photon states are based on the polarization degree of free- 
dom Q. One important example is given by states built 
on many photons [lOj and/or many qubits, and by sev- 
eral schemes of one-way optical quantum computing [11 ] . 
Hence it is of essential interest to include the use of pho- 
ton polarization in quantum circuits by fabricating in- 
tegrated polarization independent devices, i.e. able to 
efficiently guide and manipulate photons in any polariza- 
tion state. 



It has to be noticed that the above mentioned silica- 
on-silicon and UV written integrated waveguides suffer 
from intrinsic birefringence (usually reported in the order 
of 4 x 10 -4 H, [HI). In fact, these waveguides are fab- 
ricated in a doped silica multilayer structure on a silicon 
substrate and this causes material stress due to lattice 
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FIG. 2: Setup for the quantum optics experiments showing 
the source of polarization-entangled photons, the ULWBS and 
the detection system. Polarizing beam splitter (PBS), half 
(HWP) and quarter (QWP) waveplates were optionally in- 
serted in path Jca and ks to prepare different input states. 
A delay line Ax in the Ua arm enabled temporal delay varia- 
tion between the two input photons. The components shown 
in the dashed box were inserted only during the tomography 
measurement of the filtered state. C: crystal compensators, 
PC: polarization controllers. 



mismatches between the different layers. Techniques for 
reducing this stress and the induced birefringence have 
been proposed, but they pose serious difficulties in terms 
of fabrication complexity and reproducibility 
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birefringence causes polarization-mode dispersion and re- 
sults in polarization dependent behavior of the integrated 
devices, which removes indistinguishability between the 
two polarizations. Moreover, propagation in birefrin- 
gent structures can cause decoherence of large-bandwidth 
(short coherence time) photons typically generated in 
parametric down-conversion experiments. As a conse- 
quence, the techniques already employed for producing 
path-encoded quantum circuits are not appropriate for 
processing polarization-encoded qubits in integrated de- 
vices. 

In the present paper we show how to guide and ma- 
nipulate photons in any polarization state by adopting 
a recently introduced technique, based on the use of ul- 
trashort laser pulses, for direct writing of photonic struc- 
tures in a bulk glass 15, 16|. Precisely, here, for the first 
time, we demonstrate the maintenance of polarization en- 
tanglement and Bell-state analysis in an integrated sym- 
metric (50/50) beam splitter, opening the way to the use 
of polarization entanglement in integrated circuits for QI 
processes. 

Direct fabrication of buried waveguides in glass is ob- 
tained by femtosecond laser micromachining. Femtosec- 
ond infrared pulses, focused into the substrate using a 
microscope objective, induce non-linear absorption phe- 
nomena based on multiphoton and avalanche ionization. 
These processes lead to plasma formation and energy ab- 
sorption in a small region confined around the focus, 
causing a permanent and localized modification of the 
bulk material. Adjusting the processing parameters a 
smooth refractive index increase can be obtained and 
light-guiding structures are produced by translating the 



substrate with respect to the laser beam. Microscopic 
mechanisms leading to refractive index increase are com- 
plex and include densification, structural modification, 
color centers formation, thermal diffusion and accumu- 
lation. They concur in different ways depending on the 
specific material and fabrication parameter combination, 
i.e. wavelength, duration and energy of the laser pulses, 
repetition rate, objective numerical aperture (NA) and 
translation speed [17j . 

The ultrafast laser writing approach has several ad- 
vantages: i) it is a maskless technique, thus particularly 
suited for rapid prototyping of devices; ii) it can easily 
fabricate buried optical waveguides in a single step; iii) it 
can produce optical circuits with three-dimensional lay- 
outs; iv) it can provide waveguides with a circular trans- 
verse profile [18], that can support the propagation of 
Gaussian modes with any polarization state, with very 
low waveguide form birefringence. Ultrafast laser writ- 
ten (ULW) waveguides in fused silica substrate have been 
recently employed for quantum optics experiments, still 
with path-encoded qubits [l9|. However, it is known that 
ULW waveguides in fused silica are affected by material 
birefringence [20| (in particular when high refractive in- 
dex changes are required, as in the case of curved waveg- 
uides) due to the formation of self-aligned nanogratings 
in the material during the irradiation process [2l| . More- 
over fabricating waveguides in fused silica is a rather 
slow process (in the order of 10 — 100/im/s) [22]. For 
these reasons we chose to employ a borosilicate glass (EA- 
GLE2000, Corning) as substrate, where the formation of 
nanogratings has never been observed [22]. In addition, 
high repetition rate laser pulses induce isotropic ther- 
mal diffusion and melting of the material around the fo- 
cal point [23] , providing almost circular waveguide cross- 
section without the need for any shaping of the writing 
laser beam. Very low-loss waveguides are obtained with 
translation speeds as high as 1-5 cm/sec, allowing ex- 
tremely short processing times {24j|. This represents an 
advantage for the realization of complex photonic cir- 
cuits. 

At wavelengths around 800nm the waveguides support 
a single Gaussian mode of circular profile with 8/rni di- 
ameter at 1/e 2 (see the near-field intensity profile of the 
guided modes in the lower inset of Fig. [Tj), allowing a 
85% overlap integral with the measured mode of the fiber 
used (Thorlabs SM800-5.6-125) and leading to 0.7 dB 
estimated coupling losses. Measured propagation losses 
are 0.5dB/cm and using a curvature radius of 30mm ad- 
ditional bending losses are lower than 0.3dB/cm. The 
birefringence of the ULW waveguides has also been char- 
acterized {24j|, providing a value B = 7 x 10 -5 , thus 
about one order of magnitude lower than silica-on-silicon 
waveguides. 

Ultrafast laser written beam splitters (ULWBS) were 
fabricated with the directional coupler geometry, as 
shown in Fig. [TJ Straight segments and circular arcs of 
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FIG. 3: a) HOM dip with input state \HH). b) peak/dip cor- 
responding to the singlet /triplet input state, c) Fringe pattern 
obtained by rotating the HWP on mode ks- d) Fringe pat- 
tern obtained by rotating the QWP on mode kA- All curves 
represent experimental fits. 

30mm radius were employed for an overall device length 
of 24 mm. Waveguides starts with a relative distance of 
250/im and in the interaction region they get as close as 
7/im (see upper inset in Fig. [T]). This distance is the 
smallest one avoiding overlap between the two waveg- 
uides. This choice was made to minimize the sensitivity 
to fabrication imperfections and to obtain the shortest 
possible interaction length, given that future quantum 
optic devices will require several cascaded components 
integrated in the same chip. In order to optimize the 
length L of the central straight segments several direc- 
tional couplers have been fabricated varying such length 
(L = -r- 1000/im) and the corresponding splitting ratios 
were measured. L = 0/im is the shortest length yielding 
a splitting ratio of about 50% (see the ULWBS output 
modes in the lower inset in Fig. [T|) at 806nm wavelength. 
Indeed, the possibility of achieving a 50% splitting with 
no straight segments is due to the coupling between the 
modes already occurring in the curved parts of the two 
approaching/departing waveguides. The reflectivity of 
the ULWBS for the horizontal and vertical polarizations 
was measured with a tunable laser operating at 806nm. 
The measured unbalancement between the two reflec- 
tivities R H = (49.2 ± 0.2)% and R v = (58.1 ± 0.2)%, 
is attributed to a residual ellipticity in the waveguide 
profile, notwithstanding the thermal mechanism of the 
waveguide formation. Work is in progress to further op- 
timize the waveguide cross-section with the astigmatic 
beam shaping technique [Hj. 

We demonstrated the ability of the chip to preserve 
any incoming polarization state by measuring the polar- 
ization degree (G) and obtaining G > 99.8%. 



The suitability of the ULWBS to handle polarization 
encoded qubits was demonstrated by manipulating po- 
larization entangled states. The four Bell states \ip ) = 
^(\H)a\V) b ± \V) A \H) B ), \<p) = ±(\H)a\H) b ± 
I^OaI^Ob) represent an entangled basis for the four di- 
mensional Hilbert space describing the polarization of 
two photons. They can be grouped into the singlet state 
|^~), generating the antisymmetric subspace, and the 
triplet states, {|^ + )? \4>~)}-> that generate the sym- 

metric subspace, where the symmetry is referred to the 
exchange of the two photons [25[ . The beam splitter can 
be used to discriminate between the symmetric and anti- 
symmetric subspaces. Indeed if two photons in the singlet 
state \^~) impinge simultaneously on a 50/50 BS, they 
will always emerge on different outputs of the BS due to 
quantum interference. Conversely, for any state orthogo- 
nal to (thus belonging to the symmetric subspace) 
the two photons will be found in the same output mode. 

The setup adopted in the experiment is shown in Fig. 
O To observe the appearance of the bosonic coalescence 
for input symmetric states, we varied the relative delay 
between the two photons and hence their corresponding 
temporal superposition on the ULWBS [26|. We first 
tested the Hong-Ou-Mandel (HOM) [26[ effect with sep- 
arable states, by placing two PBS in the kA and ks 
modes (see Fig. [2]). We report in Figure [3^l) the co- 
incidence counts as a function of the slit delay Ax for 
two photons in the input state \HH). The measured vis- 
ibility, defined as V exp = | C °~c Q int I? where Co and d n t 
correspond respectively to the coincidence rate outside 
interference (i.e. with Ax larger than the photon co- 
herence length) and inside interference (Ax = 0). The 
measured visibility is V = 0.937 ± 0.009. We performed 
the same measurement with the input states \VV) and 
| + +) obtaining V = 0.926 ±0.012 and V = 0.954 ±0.011 
respectively. We also tested the interference with entan- 
gled states. When the photons arrive simultaneously on 
the ULWBS (Ax = in the figure) we measured for the 
triplet (singlet) a dip (peak) in the coincidence counts as 
expected (see Figure Sb))- The measured visibility are 
V\ sing i et ) = 0.930 ± 0.005 and V\ triplet) = 0.929 ± 0.005. 
We attribute the slight discrepancy observed between the 
theoretical 1 and experimental values to a partial spectral 
distinguishability between the photons on modes kA and 
fc#: this could be reduced by using narrower bandwidth 
detection filters. 

Let's now analyze the behavior of the different entan- 
gled states. The temporal delay was set at Ax = 



Since our BS has different transmission (T) and reflection (R) 
coefficients for the H and V polarization, the maximum expected 

visibilities are = 2 ThTv+Rh T*v ~ 0-989 for and V + = 
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FIG. 4: Quantum state tomography of a filtered singlet state 
on the two output mode kc and kn- Real (a) and imaginary 
(b) part of the experimental density matrix of the filtered 
state. 

and the source was tuned to generate the entangled state 
|^ + ). By inserting on mode a half waveplate (HWP) 
with the optical axis oriented at an angle 6 with respect 
to the vertical direction, the following states are gener- 
ated: — cos2#|'0 - )+sin2#|</> + ). In this case, the expected 
coincidence rates between detectors Dc and Dp after the 
beam splitter is A/o[l + V cos 40] where the expected vis- 
ibility with the given Rh and Ry can be found to be 
V = 0.973. The experimental results are shown in Fig. 
Et), yielding a visibility V = 0.962 ± 0.018. 

When the source is tuned to generate the entangled 
state \r) = -±=(\H) A \V) B - i\V) A \H) B ), by using a 
quarter waveplate (QWP) rotated at 0' on mode k Al 
the state is found to be cos 2 ^^ - ) — i sin 2 9 f \ip + ) + 

^sin2(9 / |0 i ) with If) = ±(\HH)+i\VV)). For the 
sake of simplicity, by assuming polarization independent 
reflectivity, the coincidence rates expected with a beam 
splitter with reflectivity R is A/o[l — Vteo + ^Vteo cos 4 0'\ 
where V teo = 2(1 - R)R/(2R 2 -2R+1). By taking R = 
(R H +Rv)/2 we will expect V teo = 0.987. The theoretical 
behavior was verified in the experiment. The correspond- 
ing fringe pattern, with visibility V = 0.951 ± 0.008, is 
shown in Fig. [3]i). These results demonstrate the high 
overlap between the interfering modes k A and ks and 
show that ULWBS may be used as an appropriate tool 
for the manipulation of polarization encoded qubit. 

As a final experimental characterization we adopted 
the ULWBS to carry out the projection on the singlet 
subspace. We injected into the ULWBS the separable 
state \H)a < S>\V)b and analyzed the output state when 
two photons emerge on the two modes kc and kn- the ex- 
pected state reads W~)cd- We performed the quantum 
state tomography [27] of the output state conditioned to 
the detection of the two photons in different outputs. In 
this case two standard polarization analysis setup were 
adopted after the ULWBS (see dashed box in Fig. [2j). 
The experimental density matrix, pcD, shown in Figure 
SI exhibits a low entropy (S L = 0.071 ± 0.018), a high 
concurrence (C = 0.941 ± 0.015) and a high fidelity with 
the singlet state (F = 0.929 ± 0.007). We observe that 



the present scheme achieves a-posteriori singlet compo- 
nent filtration, i.e., conditioned to the detection of one 
photon per output mode. Recently, a heralded entan- 
glement filter, based on two auxiliary photons and an 
inter ferometric scheme, has been reported by adopting a 
bulk optical scheme in Ref. [28| . 

In summary, we reported on the realization and 
quantum optical characterization of a femtosecond-laser- 
written directional coupler, acting as an integrated beam 
splitter. The experimental results demonstrate the suit- 
ability of this method to manipulate qubit encoded in 
the polarization of photon states. In order to achieve a 
complete handling of the polarization degree of freedom 
the next step will consist in the realization of integrated 
tunable waveplates and polarizing beam splitters. By 
combining these tools with integrated sources of photon 
pairs and, possibly, with integrated detectors, the real- 
ization of pocket quantum optics lab, available for optical 
quantum sensing, computing and metrology may become 
a reality in the near future. 
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amento Ateneo 2009 of Sapienza Universita di Roma and 
Finanziamento Ateneo 2009 of Politecnico di Milano. 
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DEVICE FABRICATION 

In our waveguide writing setup we employed a 
YbiKYW cavity-dumped mode- locked oscillator de- 
livering 300fs, 1/iJ pulses at 1030nm wavelength, with 
1MHz repetition rate. Laser pulses of 240nJ were fo- 
cused into the EAGLE2000 (Corning Inc.) substrate us- 
ing a 0.6 NA microscope objective. Sample translation at 
a constant writing speed of 40mm/s was enabled by high 
precision three-axes air bearing stages (Aerotech Fiber- 
Glide 3D). 

ULWBS devices were fabricated directly buried inside 
the glass substrate at a constant depth of 170/im. 

Permanent fiber coupling of the ULWBS was obtained 
by gluing two pairs of fibers at the input and at the out- 
put ports and this introduced additional coupling losses 
lower that 0.6 dB /facet. Each pair of fibers was first 
inserted into a quartz double-ferrule (Friedrich & Dim- 
mock, inc. - NJ U.S.A.), which allows to maintain a fixed 
distance of 250/im between the two fibers. 



WAVEGUIDE BIREFRINGENCE 
CHARACTERIZATION 



The birefringence of the fabricated waveguides was 
characterized by the setup shown in Fig. [TJ Six different 
polarization eigenstates (\H), \ V), |+), |— ), \L), \R)) were 
launched into the waveguide. The eigenstates were se- 
lected by rotating HWP1 and, if necessary, by adding 
also QWP1. For each of them the projections of the out- 
put state on all the same six eigenstates were measured 
(acting with waveplates HWP2 and QWP2), thus allow- 
ing the calculation of the normalized Stokes vectors for 
the output state. 

It was assumed as hypothesis that the waveguide is anal- 
ogous to a uniaxial birefringent material with the optical 
axis tilted by an angle 0, which introduces a dephasing 
5 between the two polarizations corresponding to the or- 
dinary and extraordinary refractive index. As a conse- 
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FIG. 1: Setup for birefringence measurement. LD: laser diode 
emitting at 806nm wavelength and with bandwidth A A = 2 
nm. PBS1 and PBS2: polarizing beam splitters. HWP1 and 
HWP2: half-wave plates. QWP1 and QWP2: quarter-wave 
plates. OBJ1 and OBJ2: objectives for injecting light into 
the waveguide (WG) and collecting transmitted light at the 
output. PM: power meter (PM). 



quence, a normalized Mueller matrix of the form: 



o oo 

cos 2 26 + sin 2 20cosS sin20cos20 (1 — cosS) —sin20sinS 

sin26cos20 (1 — cosS) sin 2 20 + cos 2 20cosS cos20sinS 
sin20sin5 —cos20sinS cosS 



was fitted to the experimental data in order to determine 
the transfer function from the input to the output Stokes 
vectors. S and were the free parameters in the fitting 
process. 

However, a single measurement of 5 cannot be used 
for direct determination of birefringence, since ambigu- 
ity occurs with multiple-order birefringence. Thus, the 
sample was cut to different lengths and S and were mea- 
sured as described above for each length. This allowed to 
remove the ambiguity and yielded the actual waveguide 
birefringence. 

In the ULW waveguides used for the experiments we 
found a fast axis aligned with the TE polarization and a 
birefringence value B = n e ff^M — n e ff,TE = 7 x 10 -5 . 



SOURCE 

Polarization entangled photon pairs were generated via 
spontaneous parametric down conversion in a 1.5mm f3- 
barium borate (BBO) crystal cut for type-II non-collinear 
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phase matching pumped by a CW laser diode with 
power P = 50mW and low coherence time (r p < lps). 
Degenerate photons at wavelength A = 806nm were de- 
tected within a spectral bandwidth AA = 6nm, as de- 
termined by the interference filters. Waveplates and po- 
larizing beam splitter (PBS) were optionally inserted in 
path kA and ks to prepare different input states. The 
setup shown in the dashed box have been only inserted 
in the quantum state tomography of the filtered state. 
Suitable BBO crystal compensators (C) and polarization 
controllers (PC) were respectively used on each photon 
path to compensate temporal walk-off and ensure polar- 
ization maintenance of the photons on the ULWBS. We 
could measure a coincidence rate C source = 3kHz when 
photon radiation was coupled into two single mode fibers 
and directly connected to detectors. A delay line Ax in 
the kA arm enabled temporal delay variation between the 
two input photons. 

DATA ANALYSIS 

The data shown in Figure 4 were obtained by adding 
the measurement device shown in the dashed box placed 



on the right of Figure 2, as explained in the manuscript. 
The additional optical component (HWP, QWP and 
PBS) and a more critical polarization compensation are 
mainly responsible of the measured fidelity of the recon- 
structed state. Let's try to evaluate all the imperfection 
contribution. We may attribute the 7% missing fidelity 
to the following: ~ 3% are due to the imperfection in the 
source and integrated chip (a fidelity of 97% would be 
compatible with the average visibility reported in Figure 
3), ~ 2% of missing fidelity is due to a non perfect po- 
larization compensation (in fact, the phase between the 
states \HV) and \VH) of the reconstructed state has a 
residual imaginary component as shown in Figure 4) and 
~ 2% may be attributed to the tomographic measure- 
ment setup. 
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